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Phase transitions can modify quantum behaviour on mesoscopic scales and give access to new 
and unusual quantum dynamics. Here we investigate the superfluid properties of a rotating two- 
component Bose-Einstein condensate as a function of changes in the interaction energy and in 
particular through the phase transition from miscibility to immiscibility. We show that the break¬ 
ing of one of the hallmarks of superfluid flow, namely the quantisation condition on circulation, is 
continuous throughout an azimuthal phase separation process and displays intriguing density dy¬ 
namics. We find that the resulting currents are stable for long times and possess a phase boundary 
that exhibits classical solid body rotation, despite the quantum nature of superfluid flow. To sup¬ 
port this co-existence of classical and quantum behaviour the system develops a unique velocity flow 
profile, which includes unusual radial flow in regions near the phase boundary. 


Phase transitions in quantum systems can have a dra¬ 
matic impact on the quantum mechanical behaviour on 
mesoscopic scales. Superfluidity in Bose-condensed gases 
is a mesoscopic manifestation of quantum mechanical ef¬ 
fects and one of its hallmarks is the existence of quantised 
flow around phase singularities as a response to external 
rotation mM- However, as the quantisation condition 
arises from the requirement of the single-valuedness of 
the wavefunction, an interesting, and less well investi¬ 
gated, generalization appears in superfluids composed of 
several components. In these systems, due to the in¬ 
terplay of intra- and inter-component interactions, the 
spinor order parameter can undergo a phase transition 
that modifies the global symmetry of the system. As the 
quantisation condition applies to each component inde¬ 
pendently, the path along which circulation is determined 
consequently depends on the presence of the other com¬ 
ponent. This has proven to be particularly striking in 
toroidally trapped binary mixtures of BECs, where im¬ 
miscibility can drive a transition to azimuthal phase sepa¬ 
ration, breaking the requirement of quantised circulation 
around the toroid [5] . Below we show that this transition 
is continuous and leads to a phase boundary, which ro¬ 
tates as a classical solid body. While this might seem at 
first to be incompatible with the quantum nature of su¬ 
perfluid flow, this co-existence can be explained through 
the presence of a radial flow. 

In superfluids the circulation around a closed path p 
is quantised according to v • dr = n 27 rh/m . Here n is 
an integer winding number, m the atomic mass, h the re¬ 
duced Planck constant, and the superfluid velocity field, 
V = USJOjm, is completely determined by the gradient of 
the condensate phase, 0. This implies the velocity field 
of a vortex has a tangential 1/r velocity profile, in con¬ 
trast to classical rigid-body rotation, where v = B x r. 
The creation of vortices is a response to external rotation 
and depends, in particular, on the confining geometry. 
While in simply connected trapping potentials vortices 
with higher winding numbers are unstable [6], multiply 
connected geometries are known to support persistent 


currents with large angular momentum. 

A simple multiply connected potential can be realised 
by a toroidal trap, which has recently been the subject of 
intense experimental interest [ 318 ]. For single component 
condensates superflows have been shown to exist for up 
to 40 s [9], however the superflow in toroidally trapped 
miscible two-component condensates has only been ob¬ 
served on much shorter timescales m- 

In this work we first study how the quantisation of 
circulation breaks down in the transition region between 
miscibility and phase separation in a rotating, toroidally 
trapped two-component condensate. While deep in the 
phase separation regime angular momentum scales lin¬ 
early with rotation frequency, close to the phase sepa¬ 
ration point an oscillatory behaviour is found, which is 
accompanied by significant changes in the order param¬ 
eter. At the same time, the phase profile, which drives 
the superfluid flow, adjusts in a way that allows quantum 
and classical behaviour to coexist. 

The system is modelled using two coupled Gross- 
Pitaevskii (GP) equations, which aptly describe a two- 
component Bose-Einstein condensate in the limit of zero 
temperature. Each component is assumed to be a differ¬ 
ent hyperfine state (j = 1, 2) of the same atomic species, 
mi = m2 = m, uniiii] and to consist of the same num¬ 
bers of atoms, Ni = N 2 = N. The two-dimensional cou¬ 
pled GP equations for the wave functions fjj under rota¬ 
tion around the 2 ;-axis with rotation frequency Cl = Ulz, 
are then given by 

+ E ■ rj V-i ■ 

( 1 ) 

In order to allow for stable systems with high angular mo¬ 
mentum, we assume that the atoms are trapped in a har¬ 
monic ring-shaped potential of the form Vj = \muJ^{r — 
ro)^, where vq is the toroidal radius, = x‘^ and ujr 
is the radial trapping frequency, identical for both com¬ 
ponents. Toroidal trapping potentials have recently been 
realised in several experiments by using, for example, a 
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time-averaged harmonic potential with a Gaussian laser 
beam through its centre [ 7 ] or all optical traps made by 
applying a red-detuned Laguerre-Gauss mode of a laser 
beam laiiQiiii]. Alternatively, ring-shaped trapping po¬ 
tentials which vanish asymptotically can be produced in 
the evanescent field of a nanofiber [laini. 

The coupling constants, gij = {maz\ de¬ 

scribe atom-atom interactions in terms of the three- 
dimensional scattering length and the characteris¬ 
tic harmonic oscillator length in the z direction, — 
^Jh/muJz^ For simplicity we choose the s-wave scatter¬ 
ing lengths within each component to be equal, that is 
ail = <^22 = and for both species to experience the 
same out-of-plane trapping frequencies, ujz- The strength 
of atom-atom interactions between the two components, 
^12, will be varied to induce the phase transition. 

Homogeneous two-component condensates are misci¬ 
ble for values of g‘l2 < gii922 and immiscible or phase- 
separated when g‘l2 > 911922 [HHH]- For trapped con¬ 
densates these values are slightly shifted due to the in¬ 
homogeneous density profile imilg] and the density dis¬ 
tribution in the phase separated regime is determined by 
the shape of the external trapping potential. In narrow 
ring traps (and when 9ii = 922)^ azimuthal phase sep¬ 
aration is favoured (see for instance O [ 20 l HT]), while 
in wider toroidal traps concentric ring configurations can 
occur (see for instance [ 20 l[ 2 T]). In the following we nu¬ 
merically solve the coupled GP equations by applying a 
pseudo-spectral second order Strang method with sym¬ 
metric three-operator splitting [ 22 ] . 

When the atom-atom interaction is chosen so that 
the two superfiuid components are miscible, each com¬ 
ponent is multiply connected and circulation around 
the toroid is quantised. This implies that for each 
component the average angular momentum per parti¬ 
cle, (Lz) = dx?/;* {ydjdx — xd/dy)'ipj^ is also quan¬ 

tised. In stark contrast, azimuthally phase-separated 
states break the multiply connected nature of each con¬ 
densate component around the toroid and it was recently 
shown that they can therefore rotate with arbitrary cir¬ 
culation and angular momentum of any value [ 5 ] . 

As experimentally realistic toroidal condensates are in¬ 
herently of finite size, the phase-transition takes the form 
of a continuous cross-over and in Fig. we show how the 
breakdown of the quantisation condition across this tran¬ 
sition for the above system develops. To do so, we cal¬ 
culate the angular momentum of the stationary state for 
three values of inter-atom interaction, ^12, selected such 
that the system is either fully miscible, fully immiscible 
or in the transition region between these two domains. As 
expected, for gi2 well in the miscible regime, quantisa¬ 
tion of angular momentum in each component is observed 
(blue curve). When gi2 is chosen so that the conden¬ 
sate exhibits clear azimuthal phase-separation, angular 
momentum can be seen to scale linearly with 91 (black 
curve). In the intermediate regime (^12 — ^), however. 



FIG. 1 : Phase diagram mapping the transition from quantised 
steps of angular momentum to linear scaling of {Lz)lh with 
by varying gi 2 from the miscible to immiscible regime. For 
gi 2 on the border of miscibility to immiscibility, an oscillatory 
behaviour, damped with increasing H is observed. Simulation 
parameters: = 27r x 30000 Hz, ro = 0.3 /xm with Ng = 

1.1564 X 10 “^^ (Js)2kg-h 



1 e = 27rx 700Hz 2 e = 27rxl000Hz 3 e = 27rx 1400Hz 


FIG. 2: Condensate density profiles in the transition between 
the miscible and immiscible regime. Simulation parameters: 
gi2 = 1 - 01 ^, cjr = 27r X 30000 Hz, ro = 0.3 /xm and Ng = 
1.1564 X 10 “^^ (Js)2kg-h 


an interesting damped oscillatory dependence of angular 
momentum with rotation frequency is found. The damp¬ 
ing arises as a result of the increased rotation effectively 
cancelling the harmonic trapping potential, which shifts 
the value of the critical interspecies interaction strength 
at the phase separation point towards the free space re¬ 
sult {gi2 = ^)- For larger ft the chosen interspecies in¬ 
teraction strength gi2 = l.Olg therefore moves further 
into the phase separated regime and the curve becomes 
more linear. The density distributions corresponding to 
the miscible and immiscible regimes are also displayed in 
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Fig. [2 

On the border between miscibility and immiscibility 
(red curve in Fig. [^, the density distribution changes 
as a function of the rotation frequency and three ex¬ 
amples corresponding to different angular momenta are 
shown in Fig. This behaviour can be understood by 
realising that the rotational energy acquired by the con¬ 
densate is dependent on the frequency of the externally 
imposed rotation, Q. If the system is in a phase-mixed 
state at a rotation frequency that allows for an integer 
winding number, it can acquire a certain amount of ro¬ 
tational energy with increasing rotation frequency before 
it is energetically more favourable to phase separate and 
adjust the amount of angular momentum. This leads to 
the observed cycling through mixed and phase separated 
density distributions for condensates close to the phase 
boundary. 

While in the mixing regime the well known 1 /r veloc¬ 
ity profile characteristic of rotation with quantised angu¬ 
lar momentum is exhibited, it is easy to see that in the 
phase separated regime, where fractional winding num¬ 
bers appear, this needs to be modified. In fact, if each 
superfiuid component demonstrated a perfect vortex-like 
velocity profile everywhere, the phase boundaries would 
shear over time, as the atoms closer to the center of the 
potential move faster than those at the outer radial edges. 
This would lead to an increase in the interaction energy 
and consequently unstable rotation. Instead, to ensure 
that the phase boundary is always as short as possible, 
i.e. along the radial direction, the system reacts by modi¬ 
fying the velocity profile away from purely azimuthal flow 
(see Fig. [^. 

To understand the flow profile in the phase separated 
case, we decompose the velocity field into its radial and 
azimuthal velocity contributions. These correspond to 
Vr = cos((/ 9 )i;a^ -hsin((/ 9 )i;y and = — sin((/:?)r’a, + cos((/ 9 )r’y 
and are shown in Fig. Two regions where the ve¬ 
locity field exhibits distinctly unique behaviour can be 
clearly identified. In the bulk of each component, the 
flow displays the characteristic tangential superfluid ve¬ 
locity profile of the form v oc n/re^^, with e^p a unit 
vector in the direction of the azimuthal angle ip. In con¬ 
trast, in the vicinity of the phase boundary, the velocity 
field departs from a purely azimuthal profile and a radial 
flow develops. This is consistent with the fact that each 
component only reacts to the presence of the other over 
the scale of the healing length. 

The appearance of the radial flow can be understood 
by realising that, in order to have minimal length, the 
phase boundary needs to move as a classically rotating 
object. This leads to a rotation velocity proportional to 
r, meaning that at larger radii the boundary has a larger 
velocity than at smaller radii. As this is in contrast to the 
1 /r velocity profile of the superfluid in the bulk region, 
the radial flow correctly re-distributes the atoms between 
the faster flow at smaller radii and the slower flow at 
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FIG. 3: Phase profile within the condensate, with overlaid 
lines of constant phase (pink). Upper row shows the two com¬ 
ponents in the miscible regime {gi 2 = 0.95^), where the purely 
azimuthal flow is confirmed by the fact that all lines meet at a 
single point. The lower row shows the phase separated regime 
{gi 2 — 1-6^) and the presence of a radial flow component is 
indicated by the absence of a single crossing point. Simu¬ 
lation parameters: Q = 27r x 1910 Hz, oir = 27r x 8000 Hz, 
ro = 0.75//m, Ng = 1.1564 x 10“^^ (Js)^kg“^. 



- 0.8 0.2 - 0.2 0.8 
x(//m) 



- 0.8 0.2 - 0.2 0.8 
x(/xm) 


FIG. 4: Decomposition of the condensate velocity into (a) az¬ 
imuthal and (b) radial components plotted within the conden¬ 
sate edge (defined as 6% of the maximum density Sim¬ 

ulation parameters: Q = 27r x 1273 Hz, CL;r = 27r x 10000 Hz, 
ro = 0.5/xm, Ng = 1.1564 x 10“^^ (Js)^kg“^ and gi 2 = ^-‘^g- 


larger radii. Behind the boundary, the radial flow leads 
to the movement of atoms from smaller to larger radii, 
while at the same time reducing their azimuthal velocity, 
while in front of the phase boundary the opposite process 
takes place, with atoms flowing from larger to smaller 
radius (see Fig. [^. 

The stable, classical solid body rotation of the phase 
boundary can be observed when the external rotation 
is switched off and the system is able to evolve without 
constraint (see supplementary video 1 ). As expected, the 
condensates rotate with an effective rotation frequency. 




















4 


^eff = ('^(^7^)5 which corresponds to the frequency of ex¬ 
ternally imposed rotation, that is ^ 1 . 

The above system is therefore an intriguing example 
where classical behaviour is displayed on a mesoscopic 
scale, despite the dynamics of the constituents being fully 
quantum mechanical. Similar behaviour can be found in 
the rotation of Abrikosov vortex lattices [ 23 ] . The emer¬ 
gence of large-scale classical behaviour in systems com¬ 
posed of quantum vortices also occurs in quantum tur¬ 
bulence, which displays classical Kolmogorov scaling on 
length scales larger than the average inter-vortex spacing 
[ 24 H 28 ]. 

In conclusion, we have studied the transition between 
miscibility and phase separation in rotating toroidally 
trapped two-component condensates. In the phase sep¬ 
arated regime the requirement of quantisation of circu¬ 
lation is broken and azimuthally phase separated super¬ 
fluids can rotate with arbitrary circulation. However, 
to minimise the energy of the system, the phase bound¬ 
ary has to always be aligned in the radial direction and 
therefore rotates as an effective solid body within the 
two-component flow. To resolve the dichotomy between 
this solid-body rotation, which has a velocity profile pro¬ 
portional to r, with the superfluid vortex profile in the 
bulk of the components, which is proportional to 1 /r, the 
system develops an unusual flow pattern involving radial 
components. This novel demonstration of the coexistence 
of classical and quantum behaviours can be observed in 
current state of the art cold-atom experiments. 
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